Relaxation properties of the quantum kinetics of carrier-LO-phonon interaction in 

quantum wells and quantum dots 
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The time evolution of optically excited carriers in semiconductor quantum wells and quantum 
dots is analyzed for their interaction with LO-phonons. Both the full two-time Green's function 
formalism and the one-time approximation provided by the generalized Kadanoff-Baym ansatz are 
considered, in order to compare their description of relaxation processes. It is shown that the two- 
time quantum kinetics leads to thermalization in all the examined cases, which is not the case for 
the one-time approach in the intermediate-coupling regime, even though it provides convergence to 
a steady state. The thermalization criterion used is the Kubo-Martin-Schwinger condition. 
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I. INTRODUCTION 

The Boltzmann equation is of central importance in 
many fields of physics and, since its original formula- 
tion in the theory of gases it has received a whole range 
of extensions to other domains like plasma physics, nu- 
clear physics, or semiconductors. In these fields, Boltz- 
mann scattering integrals are extensively used to model 
relaxation and thermalization processes. Adapted ver- 
sions of the H-theorem ensure that, indeed, the equa- 
tions describe the steady evolution of the system to- 
wards the proper thermal equilibrium. Early attempts 
to derive this irreversible behavior from the quantum- 
mechanical evolution have shown^ that the range of va- 
lidity of Boltzmann- like equations correspond to the low- 
coupling, slowly-varying, long-time regime. 

In more recent years with the experimental possibility 
to produce and control transport and optical phenom- 
ena at ultra-short timescales, quantum-kinetic theorieaS*^ 
have been devised in order to describe rapid processes in 
which coherence is still present, together with the onset 
of dephasing and relaxation. This means that the kinet- 
ics has to describe not only real-number quantities like 
occupation probabilities, but also complex, off-diagonal 
density-matrix elements, and their interference effects. 

Not only fast dynamics, but also the necessity to ex- 
tend the theory beyond the weak-interaction limit has 
prompted the development of quantum kinetics. A typ- 
ical example is provided by the interaction of carriers 
with LO-phonons in semiconductor quantum dots, where 
a phonon bottleneck is predicted by the Boltzmann re- 
sult (see Ref. 4 and references therein), in contrast to 
the quantum-kinetic treatment of quantum-dot polarons 
in the strong-coupling regime^ and many experimental 
findings. 

The quantum-kinetic theory using non-equilibrium 
Green's functions (GF) is one of the basic tools in this 
field. Its central object is the one-particle, two-time GF, 
for which closed equations are provided. Unfortunately, 
the large numerical effort needed for solving these equa- 



tions has limited previous applications of the two-time 
formalism to the early-time regime. The method has 
been used to describe the ultrafast optical excitation 
of semiconductors where the interaction of carriers with 
LO-phonons^, the Coulomb interaction of carriers?"'^, and 
their combined influence^ have been studied. Calcula- 
tions based on the two-time formalism also have been 
applied in plasma physics^S and for nuclear matter™. 

Since the physically relevant information (e.g. popula- 
tion and polarization dynamics) is contained in the one- 
time GF (the two-time GF at equal times), it is clear 
that a closed equation for this quantity would greatly 
simplify the procedure. This explains the huge popular- 
ity of the generalized Kadanoff-Baym ansatz (GKBA)i2i, 
an approximation which expresses the two-time GF in 
terms of its one-time component. 

The GKBA has been extensively used in the past for 
a description of non-Markovian contributions to ultra- 
fast relaxation and dephasing processes. Signatures of 
non-Markovian effects have been investigated for the in- 
teraction of carriers with LO-phonons Furthermore, 
the built-up of screening has been studied on the basis of 
a quantum-kinetic description using the GKBAi^ti^ and 
included in scattering calculationsiLi^. Results of the 
one- and two-time formulation have been compared for 
early times addressing the carrier-carrier scatteringiS as 
well as the interaction of carriers with LO-phononsiS. 

Boltzmann-like kinetic equations are obtained from the 
one-time theory based on the GKBA by further approx- 
imations: memory effects are neglected (Markov limit) 
and free particle energies are used (low coupling limit). 
One encounters therefore a situation in which only af- 
ter taking two major approximation steps, one reaches a 
kinetic theory for which the physically expected correct 
relaxation behavior can be proven analytically. To our 
knowledge, there is no attempt in the literature to explore 
systematically the relaxation properties of either the two- 
time formalism or its one-time approximation, despite 
their wide applications and the obvious fundamental im- 
portance of the problem. For example, the interest in 
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laser devices based on quantum wells^S and quantum 
dota^ requires a good understanding of the long-time 
behavior of the carriers in their evolution to equilibrium. 
Furthermore, the importance of non-Markovian effects in 
the quantum-kinetic treatment of optical gain spectra for 
quantum-dot lasers has been discussed recentlyi22iSi 

In this paper, the relaxation properties in the long- 
time limit are compared for the one-time and two-time 
quantum kinetics. As a test case, we consider the in- 
teraction of carriers with LO-phonons in semiconductor 
nanostructures, which is the dominant relaxation mecha- 
nism for low carrier densities and elevated temperatures. 
We study the optical excitation of quantum wells and 
quantum dots with short laser pulses and calculate the 
dephasing of the coherent polarization together with the 
relaxation and thermalization of the excited carrier pop- 
ulations. The equilibrium state of the interacting system 
is defined by the Kubo-Martin-Schwinger condition. We 
investigate if and under which conditions this equilib- 
rium state is reached in the time-dependent solution of 
the quantum kinetic models. This provides a unique way 
to address the range of validity of the involved approxi- 
mations. 



II. RELAXATION PROPERTIES OF THE 
BOLTZMANN EQUATION 

The Markovian limit of the kinetics, as described by 
the Boltzmann equation, is a good example to start with, 
because its relaxation properties are well understood and 
rigorously proven. To be specific, we consider the Hamil- 
tonian for the interacting system of carriers and phonons, 

^c-ph = ^ eialai + ^ hugb^bq 

i q 

+ + , (1) 

where i,j are indices for the carrier states and the mo- 
mentum q is the phononic quantum number. The cor- 
responding creation and annihilation operators for car- 
riers and phonons are given by a|,ai and b^g,bq, respec- 
tively. The Boltzmann equation for the time evolution of 
the average occupation number (population distribution) 
fi — (a|ai) has the form 

i 

with the transition rates given by Fermi's golden rule 

W.,j^^Y.\^UAq)\' (3) 
X {NgS{e, - ej - hujg) + {Ng + l)6{ei - + hug)} . 

For a phonon bath in thermal equilibrium, Ng is a Bose- 
Einstein distribution with the lattice temperature, and 



the 5-functions ensure the strict energy conservation in 
the J — >■ i transition process assisted by either the ab- 
sorption or the emission of a phonon. 

The following properties of Eq. (O can be analytically 
proven: (i) the total number of carriers fi is con- 
served, (ii) positivity is preserved, i.e., if at t = one 
has fi^O then this remains true at any later time, (iii) 
the Fermi distribution fi = [e~''('»~'^) -I- is a steady- 
state solution of Eq. (|2J) and (iv) this steady state is the 
large time limit of the solution fi{t) for any positive ini- 
tial condition provided a certain connectivity property 
holds. This property is fulfilled if any state of the carrier 
system can be reached from any other state through a 
chain of transitions having non-zero rates. The temper- 
ature of the stationary Fermi distribution is the lattice 
temperature, and the chemical potential is fixed by the 
total number of carriers. If the set of carrier states is not 
connected in the above sense, any connected component 
behaves like a separate fluid and reaches equilibrium with 
its own chemical potential. 

As satisfying as this picture looks, several problems 
arise here. The carrier-phonon interaction is essential as 
a relaxation mechanism but the carrier energies them- 
selves are taken as if unaffected by it. Both in the energy 
conserving 5-functions and in the final Fermi distribution 
these energies appear as unperturbed. This corresponds 
to a low-coupling regime, which may not be valid in prac- 
tical situations. Even in weakly polar semiconductors like 
GaAs, the confined nature of the states in quantum wells 
(QWs) and even more so in quantum dots (QDs), gives 
rise to an enhanced effective interaction24. For higher 
coupling constants one expects departures from the sim- 
ple picture discussed above. Moreover, in the case of a 
strong coupling and with the inclusion of memory effects, 
neglected in the Markovian limit, the energy conservation 
is not expected to hold. Finally, and specifically for LO 
phonons, their dispersionless spectrum, associated with 
strict energy conservation turns the system into a discon- 
nected one. Indeed, each carrier can move only up and 
down a ladder with fixed steps of size fnoLo but cannot 
jump on states outside this ladder. A phonon bottleneck 
effect in QDs was predicted on these groundsi^ 



III. STATEMENT OF THE PROBLEM 

It is clear that in most practical cases one has to turn to 
quantum-kinetic treatments in which both energy renor- 
malizations and memory effects are considered. Such for- 
malisms are provided by the two-time Green's function 
kinetics or by one-time approximations to it. In view of 
the discussion of the previous section, the following ques- 
tions, regarding the relaxation properties of the quantum 
kinetics, are in order: (i) Is the particle number con- 
served? (ii) Is positivity conserved? (iii) Is the system 
evolving to a steady state? (iv) If yes, is this steady state 
a thermal equilibrium one? In what sense? 

To our knowledge, with the exception of the first ques- 
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tion, which can be easily answered afRrmatively, there 
is no definite and proven answer available in the litera- 
ture. The aim of the present paper is to investigate how 
numerical solutions of the quantum-kinetic equations for 
realistic situations behave in the discussed respects. For 
this purpose, we compare the results of the two-time and 
the one-time approach. 



and the self-consistent field, 



(6) 

The collision term in Eq. Q has different expressions for 
and G<, 



IV. TWO-TIME QUANTUM KINETICS 

In this section we specify the Hamiltonian, Eq. 
for the case of a homogeneous two-band semiconductor, 
where carriers interact with LO-phonons via the Frohlich 
coupling, 

fe,A q 

+ 4+q,Aafc.A(foq + hl^) . (4) 

fc,q, A 

The carrier quantum numbers are the band index A = c, w 
and the 3D- (for the bulk case) or 2D- (for QWs) momen- 
tum k. The coupling is defined by ^ a/q^ for the 3D 
case, or by ~ aF{q)/q for the quasi-2D case, with the 
form factor F{q) related to the QW confinement func- 
tion and the Frohlich coupling constant a. Additional 
terms to this Hamiltonian describe the optical excita- 
tion and the Coulomb interaction in the usual wayi^ We 
consider only sufficiently low excitations so that carrier- 
carrier scattering and screening effects are negligible. 
Then the only contribution of the Coulomb interaction 
is the Hartree-Fock renormalization of the single particle 
energies and of the Rabi frequency. 

The object of the kinetic equations is the two-time GF, 
G^.'^ {ti,t2) — G^'^ {t,t — r). We use the parametriza- 
tion of the two-time plane (ti,i2) in terms of the main 
time t and relative time r. One can combine the two 
Kadanoff-Baym equations^^ which give the derivatives of 
the GF with respect to ti and t2, according to d/dt — 
d/dti + d/dt2 and d/dr ~ —d/dt2 in order to propagate 
the solution either along the time diagonal (i-equation) 
or away from it (r-equation). As two independent GFs 
we choose the lesser and the retarded ones, and limit our- 
selves to the subdiagonal halfplane r ^ 0, since supra- 
diagonal quantities can be related to subdiagonal ones 
by complex conjugation. With these options and in ma- 
trix notation with respect to band indices, the main-time 
equation reads 



d 



^f^Q-Gt^{t,t- r) = ^iit) G'^^^it,t- t) 



iR.<, 



(5) 



,H-G^{t,t-T] 



d 



dt' 



j:^it,t')G^it',t~T) 



coll Jt-T 

-G^{t,t')i:i{t',t^T) 



(7) 









r) 




f dt' 




coll 


/ — oo 



(8) 



The time arguments of the self-energies and GFs in 
Eq. JSJ are the same as in Eq. Q and are omitted for sim- 
plicity. The advanced quantities are expressible through 
retarded ones by conjugation. The self-energies are com- 
puted in the self-consistent RPA scheme and have the 
explicit expressions 

= ^hY,gl [D>{t - t')G^_^{t,t') 
+ D^{t-t')G<_^{t,t')] , 

E< (t, t')^zhJ2 9l D< [t - t')G<_^ it, t') , (9) 

with the equilibrium phonon propagator 

£)|(i) = -1 Nq e^''"^"' + (1 + Ng) e^''"^"^ . (10) 
For practical calculations, we use dispersionless phonons 

— ^LO- 

The above set of equations has to be supplemented by 
specifying the initial conditions. For all the times prior 
to the arrival of the optical pulse, the system consists of 
the electron-hole vacuum in the presence of the phonon 
bath. This is an equilibrium situation, characterized by 
diagonal GFs, which depend only on the relative time. 
More precisely, one has 



G^{t,t-T) 



G<{t,t-T) = 



9?{r) 
9^[r) 







(11) 



where the instantaneous self-energy contains the external 



where is the retarded GF of the unexcited system. 
The calculation of this GF is an equilibrium problem. 
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namely that of the Frohlich polaron. The polaronic GF 
is the solution of the r-equation, 

{^h^ - el^ alxir) = £ dr' E^;,(r - r').9f;,(r') , 

(12) 

in which, to be consistent with the i-evolution described 
above, the RPA self-energy is again used. The vacuum 
GF of Eq. ((TT|l is not only the starting value for the GF 
in the main-time evolution but it also appears in the in- 
tegrals over the past which require GF values before the 
arrival of the optical pulse. Moreover, the presence of 
the polaronic GF brings into the picture the complexi- 
ties of the spectral features of the polaron, with energy 
renormalization and phonon satellites. Finally, the decay 
of the polaronic GF introduces a natural memory depth 
into the problem. An example is seen in Fig.^where, due 
to a rather strong coupling constant and a high temper- 
ature, the decay with the relative time r is rapid. This 
allows to cut the infinite time integrals of Eq. (jSJ at a 
certain distance away from the diagonal. 

In Fig. ^ the momentum argument is replaced by the 
unrenormalized electron energy E = h^k'^/2m*. This 
change of variable is allowed by the fact that the momen- 
tum dependence of the polaronic GF is isotropic. The 
same energy argument is used in the subsequent figures, 
with the exception of Fig. |21 where the reduced mass is 
employed {E = h^k'^ /2m*) as being more appropriate 
in the case of the polarisation. The choice of an energy 
variable facilitates the comparison with other energies ap- 
pearing in the theory. For instance, in Fig.^a somewhat 
slower decay is seen at low energies and is a trace of the 
phonon threshold {Hujlo = 21 meV). Also, in Figs. I2I3I 
below, the functions have a peak around the detuning 
(120meV). 




FIG. 1: Absolute value of the polaronic retarded GF for elec- 
trons in a CdTe QW at T = 300K. £ = corresponds to the 
conduction band edge. 

With the specification of the initial conditions, the 
problem to be solved is completely defined. After obtain- 
ing the two-time GFs, the physically relevant information 
is found in the equal-time lesser GF, which contains the 
carrier populations and the polarization. This program 
is carried out for the case of a CdTe (a = 0.31) QW at 
room temperature. The excitation conditions are defined 



by a Gaussian-shaped pulse of 100 fs duration (FWHM of 
the intensity), having an excess energy of 120meV above 
the unrenormalized band gap and an area of 0.05 (in tt 
units). This gives rise to carrier densities in the order of 
10^/cm^, sufficiently low to neglect carrier-carrier scat- 
tering. 



P 




E (meV) 



FIG. 2; (Color online) Time evolution of the coherent inter- 
band polarization after optical excitation of a CdTe QW with 
a 100 fs laser pulse centered at time t = 0, using a two-time 
calculation. 




FIG. 3: (Color online) Time evolution of the electron popu- 
lation distribution for the same situation as in Fig. |21 

As seen in Figs.|21and|31 the interaction of carriers with 
LO-phonons provides an efficient dephasing and leads, 
in a sub-picosecond time interval, to a relaxation of the 
electron population into a steady-state distribution. The 
same is true for the hole population (not shown). Before 
discussing this result we compare it to the outcome of 
the one-time calculation. 



V. ONE-TIME QUANTUM KINETICS 

To obtain from the Kadanoff-Baym equations for the 
two-time GF a closed set of equations for the equal- 
time lesser GF, G'^{t,t), one can use the generalized 
Kadanoff-Baym ansatz (GKBA)^^. The ansatz reduces 
the time-offdiagonal GFs appearing in the collision terms 
of Eq. JSjl to diagonal ones with the help of the spectral 
(retarded or advanced) GFs. Therefore, the GKBA has 
to be supplemented with a choice of spectral GFs. In our 
case, it is natural to use for this purpose the polaronic 
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GF. For r > 0, this leads to the GKBA in the form 

G<{t,t^T) « ihg^ir) G<{t~T,t~T) . (13) 

The resuh of this procedure for the same system and 
using the same excitation conditions as for the two-time 
calculation is shown in Fig. ^ We find that the steady 
state obtained in this way differs appreciably from that 
of the two-time calculation. 




FIG. 4: (Color online) Time evolution of the electron popu- 
lation distributions in a CdTe QW, using the same excitation 
conditions as in Figs. HEl and a one-time calculation. 



VI. THE KMS CONDITION 

For a fermionic system in thermal equilibrium, the fol- 
lowing relationship connects the lesser and the spectral 

G<{uj)^-2i Jiuj) Im G^iuo) , 
f(u;) = ^.^ ^ , . (14) 

In thermodynamic equilibrium, the GFs depend only 
on the relative time and Eq. 114|) involves their Fourier 
transform with respect to this time. The relationship is 
known as the Kubo-Martin-Schwinger (KMS) condition 
or as the fluctuation-dissipation theorem, and leads to a 
thermal equilibrium population given by 

/.^ = -/^/HlmGf-(c.). (15) 

The two-time theory provides the excitation-dependent 
retarded GF along with the lesser one, the formalism be- 
ing a system of coupled equation for these two quantities. 
Nevetheless, in the low excitation regime used here the 
difference between the actual retarded GF and its vac- 
uum counterpart \{^) turns out to be negligible, as 
can be checked numerically. Therefore, the latter can be 
used in Eq. ifT^ without loss of accuracy. 

The thermal equilibrium distribution obtained from 
the KMS condition is the generalization of the Fermi 
function of the non-interacting case and is used as a check 
of a proper thermalization. The test of the two steady- 
state solutions against the KMS distribution function is 



seen in Fig. [S| The two-time calculation is in good agree- 
ment with the KMS curve, but the one-time evolution 
is not. It appears that the one-time kinetics produces a 
steady state with a temperature considerably exceeding 
that of the phonon bath. 




50 100 150 200 250 

E (meV) 

FIG. 5: (Color online) One- and two-time CdTe QW electron 
populations at t = 1240 fs and the KMS result. 

It is to be expected, however, that for a weaker cou- 
pling the discrepancy between the full two-time proce- 
dure and the GKBA is less severe. This is indeed the case, 
as shown in Fig.|21 where results for a GaAs {a = 0.069) 
QW are given. The wiggles seen in the two-time curve 
are traces of the phonon cascade, which are still present. 
This is due to the much longer relaxation time in low- 
coupling materials. Nevertheless the trend is clear, the 
steady-state solutions of both approaches are in good 
agreement with the KMS condition. 




50 100 150 200 250 300 350 400 

E (meV) 

FIG. 6: (Color online) Electron population at t—1600 fs for a 
GaAs QW and optical excitation with a 100 fs laser pulse at 
t — 0. Solutions of the two-time and the one-time quantum 
kinetics are compared with the KMS result. Inset: same on 
semi-logarithmic scale. 

Another important example concerns a non- 
homogeneous system. It consists of CdTe lens-shaped 
self-assembled QDs, having both for electrons and for 
holes two discrete levels below the wetting-layer (WL) 
continuum. These states are labelled s and p, according 
to their z-projection angular momentum. We consider 
an equidistant energy spacing of 2Ahu}LO between the 
WL continuum edge, the p-level and the s-level, for 
the electrons and a Q.21hwLo similar spacing for holes. 
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The formalism used is the same as for the homogeneous 
systems but with the momentum replaced by a state 
quantum number running over the discrete QD states 
and the WL continuum. This amounts to neglecting 
GF matrix elements which are off-diagonal in the state 
index, but still keeping off-diagonal terms with respect to 
the band index. This has been shown to be a reasonable 
approximation for QDsi^iS& Our calculations for this 
example include both localized QD and delocalized WL 
states. We consider a harmonic in-plane confinement 
potential for the localized states and construct orthogo- 
nalized plane waves for the delocalized states in the WL 
plane. The strong confinement in growth direction is 
described by a step-like finite-height potential. Details 
for the calculation of interaction matrix elements are 
given in Ref. i5|. In Figs. [7| and |S1 the time evolution 
of the population of electrons is shown. The system is 
pumped close to the renormalized p-shell energy with a 
100 fs laser pulse at time t — 0. Therefore the majority 
of the carriers is initially found in the p-state (which has 
a two-fold degeneracy due to the angular momentum in 
addition to the spin degeneracy). Nevertheless, efficient 
carrier relaxation takes place, even if the level spacing 
does not match the LO-phonon energy, and a steady 
state is reached. The two-time results are again in 
agreement with the KMS condition, shown by open 
circles. The one-time evolution shows a non-physical 
intermediate negative value for the WL population and 
converges to a state in strong disagreement with the 
KMS resuh. 




400 800 1200 



t (fs) 

FIG. 7: (Color online) Electron populations in the localized 
s and p states of a CdTe QD and in the extended fc = 
WL state after optical excitation with a 100 fs laser pulse 
at t = 0, as calculated using the two-time kinetics. Open 
circles represent the equilibrium values according to the KMS 
condition. 



VII. CONCLUSIONS 



The long-time behaviour of different quantum kinetic 
approaches to the problem of carrier-scattering by means 
of LO-phonons was analyzed, in order to assess their re- 
laxation properties. As a test of proper convergence to 
thermal equilibrium, the KMS condition was used. We 
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FIG. 8: (Color online) Same as Fig. |7| but using the one- 
time kinetics. Note that an identical ordinate axis is used to 
facilitate the comparison. 



considered materials with low (GaAs) and intermediate 
(CdTe) Frohlich coupling. The results can be summa- 
rized as follows: (i) In both the one-time and the two- 
time quantum kinetics steady states are reached, (ii) The 
steady state produced by the two-time approach obeys 
the KMS condition in all cases considered, (iii) The one- 
time result agrees with the KMS condition only at low 
coupling and differs considerably for larger ones. 
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